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Renormalization group analysis on hierarchical lattices is a valuable tool to understand the critical 
behavior of statistical mechanical models. In presence of quenched disorder, where position space 
renormalization group procedures are still on their way, it has the advantage to be analytically ap- 
plicable. In some model cases, though, predictions obtained with the Migdal-Kadanoff bond removal 
approach fail to qualitatively reproduce critical properties obtained in mean-field approximation or 
by numerical simulations. We focus on the critical behavior of Ising and Blume-Emery-Griffiths 
models on hierarchical lattices microscopically more similar to finite dimensional Bravais lattices 
than those belonging to the Migdal-Kadanoff class in order to understand if and how different 
behaviors are associated with hierarchization. 



I. INTRODUCTION 

Understanding the nature of the low temperature 
phase of spin-glasses, i.e., systems with quenched dis- 
ordered magnetic exchange interaction, in finite dimen- 
sional systems has turned out to be an extremely dif- 
ficult task. Since the resolution of its mean-field ap- 
proximation, valid above the upper critical dimension 
(D = 6), more than thirty years have passed without 
a final word about the possible generalization of mean- 
field properties of spin-glasses to finite dimensional cases. 
The mean-field, else called Replica Symmetry Breaking 
(RSB) theory™ involves a very interesting solution for 
the spin-glass phase and its critical properties, rich of 
physical (and mathematical) implications, and has been 
fundamental in solving very diverse problems both in 
physics and in other disciplines^^ 

Because of its complicated structure, to overcome tech- 
nical (maybe also conceptual) obstacles hindering the 
"portability" of RSB theory predictions to short-range 
systems on Bravais lattice in D < 6 is a rather big chal- 
lenge in theoretical physics. Indeed, the RSB solution is 
so complex that non-perturbative effects cannot be taken 
under control in any perturbative loop-expansion around 
the upper critical dimension and critical scaling behavior 
is yet to be understoodP^ The main hindrance is the 
lack of translational invariance in the position space for 
locally frustrated systems with quenched disordered in- 
teractions, making the techniques developed for quantum 
field theory and succes sfully exported to statistical me- 
chanical problem d 12 * 13 * (e.g., for the Ising model critical 
exponents) inapplicable. 

For what concerns Wilson's original approach in po- 
sition spacepl the extension of renormalization group 
techniques to disordered and locally frustrated systems 
is still on its way. The generalization of classic position 
space renormalization methods on Bravais lattices to dis- 
ordered interaction, such as the o nes proposed for Ising 
spin models in the seventies] 15 * 16 ! has led to controver- 
sial results. On the one hand, by means of a cumulant 



expansion approach, evi dence for a spin-glass phase is 
yielded in dimension two| 17 | 18 J lower than the lower criti- 
cal dimension D = 2.5.Q2H2DOn the other hand the renor- 
malization through block transformation on spin clusters 
does not yield any spin-glass fixed point even in dimen- 
sion threeP^ 

The only consistent results have been achieved using 
"realisable" approximations,^ namely those that are the 
exact solution of some alternative problem. So most of 
the studies have been concentrating on hierarchical lat- 
tices for which, in the ordered cases, the position space 
renormalization group flow is indeed exact (no trunca- 
tion required). The study of these systems has brought 
to very important results, cf., e.g., Refs. [22 -24 and Ref. 
[2 5) and references therein. 

In this work we will investigate the renormalization 
group analysis on spin systems with quenched disorder 
on hierarchical lattices. We shall consider both Migdal- 
Kadanoff (MK) as well as more complex hierarchical lat- 
tices. 

Our main aim is to investigate whether hierarchical 
cells locally more similar to short-range Bravais lattices 
can reproduce features of spin-glasses and ferromagnets 
so far unobserved in position space renormalization group 
studies, and then try to understand the limits and the 
advantages of this approach. As byproduct we will also 
provide new estimates for critical quantities on several 
hierarchical lattices. 

As an instance, MK lattices fail to represent short- 
range spin-glasses on Bravais lattices also in the mean- 
field approximation and are, thus, strongly limit ed in 
probing the actual nature of the spin-glass phasePSE3 
For what concerns perturbative disorder in ferromagnetic 
systems, it has been found that Nishimori conjecture does 
not provid e th e exact value of the multicritical point 
coordinate, 2824 although recently Ohzeki, Nishimori e 
Berkei^^HS! have put forward an improved conjecture for 
models defined on hierarchical lattices and they found 
that for various families of these models a noteworthy re- 
covery is obtained. Furthermore, MK lattices fail, even 
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qualitatively, in obtaining first order transitions both in 
ordered^ and in disordered^ systems. Eventually they 
do not show re-entrances observed in some spin-1 spin 
systems, e.g., the Blume-Emery-Griffiths models, both 
in mean-field theory 3 ^! and by numerical simulations in 
finite dimension.^ Our aim is to test whether and which 
of the above mentioned differences are specifically due to 
the bond moving procedure at the basis of the renor- 
malization group analysis. We will, thus, implement 
and compare the analysis of the critical behavior of well 
known statistical mechanical models with quenched field 
and bond randomness on both MK hierarchical lattices 
and more complex hierarchical lattices, such as Wheat- 
stone bridge (WBJp31 r folded-celf^l lattices. 

The paper is organized as follows: in Sec. |n]we expose 
the implementation of the position space renormalization 
group we used; in Sec. |III| we investigate Ising spin mod- 
els, presenting and comparing the estimates of critical pa- 
rameters and discussing how they comply to known sta- 
tistical mechanics criteria in presence of disorder (Nishi- 
mori conjecture, Harris criterion, ferromagnetic line in- 
version, ...). We consider first, Sec. IIIB the random 
field Ising model. Then, we study the Ising model with 
quenched bond disorder. In Sec. |III C| we analyze the 
behavior in low dimension, i.e., d < 2.5 where only para- 
magnetic and ferromagnetic phases are displayed, and, 
in Sec. HI D| in d > 3, where the spin-glass phase arise. 
In Sec. |IV| we consider the Blume-Emery-Griffiths model 
on several hierarchical lattices in dimension d > 3. In 
the latter case our analysis shows a phase diagram dis- 
playing a reentrance for strong disorder, absent on MK 
lattices,^ but present in the mean-field approx imation^ 
and numerical simulations on 3D cubic lattices. ^ 1 * 37 * 38 ! 



II. HIERARCHICAL RENORMALIZATION 

The Positional Space Renormalization Group (PSRG) 
approach, approximated on realistic Bravais lattices, be- 
come exact as it is iterated on Hierarchical Lattices 
(HL) I221M139HH] These lattices are constructed by carrying 
successive similar operations at each hierarchical level. 
E.g., at each level one replaces bonds by well-defined unit 
cells (see for example Fig. [I] for the diamond lattice or 
Fig. [2] for a MK lattice). The PSRG procedure works 
the inverse way of the lattice generation, i.e., through a 
decimation of the internal sites of a given cell, leading 
to renormalized quantities associated with the external 
sites. 

In the pure case, the PSRG analysis proceeds as known 
by finding the interactions leaving the partition function 
invariant under decimation, and obtaining the critical ex- 
ponents by the eigenvalues of the first derivative matrix 
computed on the relative fixed point 

In disordered systems the PSRG trasformation is de- 
scribed by the evolution of a probability distribution 
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FIG. 1: Costruction of diamond, else called 2-dimensional 
Wheatstone-bridge, hierarchical lattice with fractal dimension 
d = log 5/ log 2 = 2.3219... . 




FIG. 2: Necklace MK lattice. It has b - 
sion d = 3. 



2 and fractal dimen- 



rather than single values of coupling constants! 43 * 44 ! 



P'(J') 



n dj a p{j a ) 



S[J' -K({JM] 



(1) 



where J is the set of external parameters (couplings, 
fields, chemical potentials, ...),& is the the length of the 
cell in lattice spacings (i.e., the scaling factor in the dec- 
imation procedure), d the space dimension, so that b d is 
the size of the cell in number of bonds to be decimated, 
and TZ({J'}f > d) is the local recursion relation for the in- 
teractions. 

In MK lattices, because of their ID-like topology, see, 
e.g., Fig. [2j the transformation can be divided into steps 
of so-called bond-moving and decimation, each of which 
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FIG. 3: RG evolution of the probability distribution of nearest-neighbor interaction for the spin-glass Ising model on the 
Wheatstone bridge lattice of Fig. [i] (D ~ 3.58) and pool size M = 10 6 . Paramagnetic (left), Ferromagnetic (center) and 
Spin-Glasss (right) phases. 



involving only two bonds at a time. It is, thus, possi- 
ble to exactly compute the probability distribution and 
represent it with histograms, each bin of which charac- 
terized by a value of the interactions and an associated 
probability.^ 

In HLs more complicated than MK, such as the 
Wheatstone-bridge-like in Fig. [T] this factorization is no 
longer possible and we must consider the convolution of 
more than two links at a time. We are, eventually, obliged 
to proceed in a statistical way. The PSRG scheme is, 
then, accomplished by representing the probability dis- 
tribution of the couplings by a pool of M real numbers^ 
from which one can compute its associated moments, at 
each renormalization step. In the limit M — > 00 these 
moments should approach those of the exact renormal- 
ized probability distribution. The process starts by creat- 
ing a pool with M coupling constants generated accord- 
ing to the initial distribution. A PSRG iteration consists 
in M operations in which one randomly picks a set of b d 
couplings from the pool in order to generate one renor- 
malized coupling, which will populate the renormalized 
pool. Following this procedure, one creates a new pool 
of size M representing the renormalized probability dis- 
tribution. During the PSRG procedure the moments of 
the coupling distribution are of particular interest for the 
identification of the phases. 

For example, in Ising models with quenched disorder, 
by denoting with J the average of the couplings and with 
a j the mean square displacement, one obtains the Param- 
agnetic (PM), Ferromagnetic (FM), and Spin Glass (SG) 
phases, as dominated by the attractors 

J^O; aj^Q; PM ; 

J — » 00; a j — > 00 (J/<Jj — > 00) ; FM ; 



J -> 0; 0-7 -> 



SG 



In Fig. [3] the typical PSRG iterations of the couplings 
distribution in the three phases are shown for the random 
bond Ising model on the HL lattice of Fig. |4j 

In order to reduce the dependence on a particular se- 
quence of random numbers, the evolution of each distri- 



bution is analyzed over Ns different samples. This is es- 
pecially relevant when the starting pool is near a critical 
point, and random fluctuations can lead different sam- 
ples of the same distribution into different attractors. In 
this case, we will adopt the convention that a phase, is 
considered as identified if at least 80% of the Ns samples 
flow into the same attractor. This defines the error for 
the location of critical points, which can be reduced by 
increasing the value of M. 



III. MODELS WITH ISING SPINS 

The well known Ising model is defined by the Hamil- 
tonian 



f3H(s) = J 2~] SiSj + h 22 s t 



(2) 



where = ±1 and (ij) indicates a sum over nearest- 
neighbor pairs. We stress that here and in the rest of the 
paper we include the temperature in the definition of the 
couplings (reduced parameters). 

Since we, eventually, want to study the critical proper- 
ties of disordered systems, and since through the renor- 
malization group transformation original random bonds 
induce random fields (and vice-versa) already at the first 
step of renormalization, it becomes more convenient to 
start using the following Hamiltonian 



E 



Jij SiSj 



Si + Si 



h 



(3) 



In this way each link between two sities i and j is asso- 
ciated with three (possibly disordered) interactions J^, 



hij and h}j. 



Decimating the inner sites {s} of the basic cell C a b of 
the hierarchical lattice with external sites s a and s&, while 
imposing the conservation of the partition function of the 
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cell 

Zc ab = x SaSb = ^ cxp{-f3H[s a ,s b ;{s}}} , (4) 
{s}ec ab 

yields the renormalization group equations: 
1 



Jr = 7 lo S 



Xa a;_ 

+ 



2 \x |- 



(5) 



The partition sums x SaSb , also called principal Boltzmann 
factors of the cell, are the weights of the cell for fixed 
external spins s a and Sb- The sum in Eq. @ runs over 
all inner or free spins of the cell Cab- 
in the zero-temperature limit the relations become 

4J K = max[-pH(l,l,s)] + max[-f3H(-l,-l,s)] 
- max[-j8K(l,-l,s)] - max[-j0W(-l, 1, s)] 
2h R = max[-^K(l, l,s)] - max[-jSK(-l, -1, s)] 
2/i^ = max[-jSK(l,-l,s)] - max[-j9W(-l, 1, s)] 

(6) 

When the external field is missing h = h) = and 
H(s) = H(s), which = h ] R = 0. 




FIG. 4: Wheatstone bridge (WB) hierarchical lattice obtained 
by bonds moving from a cubic lattice. It has 6 = 2 and fractal 
dimension d = log 12/ log 2 « 3.585. 





A. The ordered ferromagnetic Ising model 



For an ordered ferromagnetic system (Jy = J, /iy = ft, 
and h\j — hJ) the critical exponents can be obtained from 
the eigenvalues of the first derivatives matrix 



9Jr dJ R dJ R \ 

dJ dh dh1 \ 

dh R dh R dh R 

dJ dh dhi 

dh\ dh R dhj ) 

dJ dh ~dhT ' 



(7) 



computed on the pure fixed point corresponding to the 
universality class of the ferromagnetic transition. The 
derivatives are easily obtained using 



dx. 



dJ 

s a s, 

dh 

$X s a s, 



E E StS i I eX P[-^( S »- S b, {«})] 



E E 

{s}ec ab 

E fe 

{s}£C ab 



1 exp[-pH(s a ,s b ,{s})] . 
' ' exp[-P*H(s a ,Sb,{s})] 



(8) 



FIG. 5: Left hand side: square cells of cell spacing length 
6 = 3. On the right side: corresponding folded square hier- 
archical lattices with two roots (open circles). All outcoming 
and incoming sites, pointed by the arrows, are put together 
and generate, respectively, root sites a and b. The inner sites 
of the square cell become the inner sites of the folded square 
hierarchical lattice. Note that with this construction we ob- 
tain self-dual lattices with fractal dimension two. 



In particular, if the fixed point is for h = fv = 0, 
it is easy to see that the matrix in Eq. (J7]) is diagonal 
and dh^/h^ = c, where c is the number of incoming 
(outgoing) link in the external outgoing (incoming) site. 
For example c = 4 in Fig. |4j c = 3 in Fig. [3] and c = 5 in 
Fig. [6] In this case the only relevant eigenvalues are At = 
8tT' — djJ' and A^ = dhh' , with the corresponding 
scaling exponents utm = log;, Xr,h- 

The above PSRG scheme holds when the coupling con- 
stants Jij are all equal and the external field is ordered 
and homogeneous. To perform the equivalent analysis of 
the critical behavior in disordered systems one has to im- 
plement a PSRG method based on Eq. ([I]) and look for 
fixed point probability distributions. In this latter case, 
though, a general method for finding the critical expo- 
nents - like the one provided by Eqs. Q, ^ - is hard 
to be provided. The idea is, then, to estimate the criti- 
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FIG. 6: Folded square lattice obtained as in Fig. [5j but for 
6 = 5. 



cal exponents by slightly perturbing the system from the 
unstable fixed point distribution and measure how fast 
it departs from it under successive PSRG iterations. We 
will discuss this procedure in detail in the following. 



B. Random Field Ising Model 

In this section we discuss the PSRG study of the Ran- 
dom Field Ising Model (RFIM) with bimodal and Gaus- 
sian distributed quenched external field on the simple 
necklace MK lattice of Fig. [5] with fractal dimension 
d = 3, and on the Wheatstone-Bridge (WB) hierarchical 
lattice of Fig. [4j with fractal dimension d sa 3.585. For 
the bi moda l case on the MK lattice we recover known 
resultsP^H The results on the WB lattice are, instead, 
obtained for the first time. We stress that in a previ- 
ous paper of Nobre and SalmorP^the exact PSRG trans- 
formation of the hierarchical lattice is not achieved, as 
pointed out by BerkerP^ 

The initial distribution of couplings for the RFIM reads 



P(Jij 5 hi J 1 ; ^ 



)=S(J ij -l)p(h ij )S(h\ j ) 



(9) 



where p(hij) is either a bimodal or a Gaussian distribu- 
tion: 



p(hij)-- 



1 [6 (hi j- h ) +5 (hi j- 

7^! exp {"^t}- 



(10) 



The initial distribution is an even function of h, and 
h) . This symmetry is preserved under the PSRG trans- 
formation. To maintain this symmetry in our finite sam- 
ple, we, actually, use a pool of 2M interactions: we ex- 
tract M values from the pool of renormalized interac- 
tions (Jij , hij , h\j) and M from (Jij,—hij,—h\j). This 
trick becomes important for long PSRG iterations and 
as M increases the PSRG threshold step beyond which 



it becomes necessary quickly increases. The forced sym- 
metrization is, thus, not crucial in determining critical 
properties but helps decreasing finite size effects for small 
pools. For the present computation we take pools with 
up to M = 10 6 interactions, enough to yield statistically 
stable results. 

The RFIM shows two phases identified by the behavior 
of the ratio between the average coupling J = (Jij) and 
the standard deviation of the fields h\ = (hfj) along the 
PSRG flow: 

ho/J^>-oo high temperature phase; 

ho/J—>0 low temperature phase; 

(and similarly for h' field). At the critical point both J 
and ho (as well as <r h t) flow to infinity, confirming that 
the critical behavior of RFIM is controlled by a zero- 
temperature fix point couplings distribution. 

For the MK lattice of Fig. [2] the critical point is ob- 
tained for the bimodal case at ho = 0.942(4), with a 
fix point ratio h /J — 0.876(3), whereas for the Gaus- 
sian case at ho = 0.934(5), and fix point ratio ho/J = 
0.869(4). 

For the WB lattice shown in Fig. [4] the critical point 
is obtained in the bimodal case for ho = 0.460(3), with 
a fix point ratio ho/J — 0.443(5), whereas for the Gaus- 
sian case at ho — 0.456(3), and fix point ratio ho/J = 
0.448(4). 

In Fig. [7] we show the projections of the critical fix 
point probability distribution for the bimodal case on 
the WB hierarchical lattice. In the Gaussian case they 
are qualitatively the same. 

At the zero-temperature fix point there are three inde- 
pendent exponents z, y^ and yr = which can be 
obtained following the procedure of Ref. |47j for both the 
bimodal and Gaussian case. 

Distribution growth scaling. The exponent z de- 

scribes the deviation of the couplings probability distri- 
bution from the unstable fix point distribution under the 
PSRG transform. In order to estimate z, once the PSRG 
flux gets close to the fix point distribution, we fix the 
ratio ho/J at the critical value in the subsequent PSRG 
iterations by shifting at each step all the couplings {Jij} 
in the pool towards the ideal critical value. In particular, 
we choose to change them of the 80% of the difference 
with the critical value. In our computation each shift 
turns out to be very small, of the order of 0.1%, though 
essential because of the unstable nature of the fix point. 

We, then, evaluate the rescaling factor A = h' /ho = 
J/J' at each PSRG step and then take its average. The 
exponent is then estimated as 



z = logfe A , 



(11) 



where in our case the overbar denotes the average over 
10 PSRG steps, obtaining: on the MK cell, cf. Fig. 
[2] z = 1.491(3) for the bimodal and z = 1.486(3) for 
the Gaussian distribution; on the WB cell, cf. Fig. |3J 
z = 1.788(2) for the bimodal and z = 1.787(2) for the 
Gaussian distribution. We note that in all the cases it is 
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z < d/2 2 that is the upper bound provided by Berker and 
McKayP 

External field scaling. The exponent yh describes 
the rescaling of an infinitesimal homogeneous field and 
can be obtained by averaging the relations in Eq. Q, 
|8]) over the fix point distribution 

»-■*($)• <12 » 

For MK we obtain = 2.991(1) for the bimodal and 
yh = 2.990(1) for the Gaussian distribution. For WB we 
obtain yh = 3.575(1) for the bimodal and yh = 3.576(1) 
for the Gaussian distribution. In all cases the value is 
smaller than the fractal dimension d, which is 3 for MK 
and ~ 3.585 for WB, implying that the magnetization is 
continuous at the transition. 

Correlation length scaling. In order to estimate 

the exponent v, we first reach a pool of renormalized 
couplings satisfactorily representing the fix point distri- 
bution. Next we take a copy of the pool and generate a 
slightly perturbed couplings probability distribution by 
shifting every coupling {J^ } of the replicated pool by a 
small amount 5 = 10 -4 • J. The two pools, the original 
and the perturbed one, are then simultaneously renor- 
malized. To reduce statistical fluctuations the couplings 
in the pool representing the fix point probability distri- 
bution are shifted, after each PSRG step, to keep the 
distribution close to the unstable fix point, in a manner 
similar to that used for the estimation of the exponent 
z. Note that in this way the only role of the shift 5 is 
to accelerate and make explicit the departure from the 
fixed point. 

By defining t n as the difference between the value of 
the ratio ho/ J in the two pools after n PSRG iterations, 
the correlation length exponent is estimated as 



where the overbar denotes the average over the PSRG 
iterations n. Note that the argument of the logarithm is 
always positive, because leaving the fix point the second 



copy variance can either shrink or increase in its flux, but 
it does not oscillates between different PSRG steps. 

Typically we have n ~ 3 -j- 9, for which the perturbed 
pool is not too far from the unstable fix point, and the 
ratio ho/ J is stable. The result is independent of M and 
it is quite stable over independent PSRG evolutions, at 
least for M > 10 5 . 

We obtain \/v = 0.45 ± 0.23 for the bimodal distri- 
bution and \/v = 0.44 ± 0.20 for the Gaussian distribu- 
tion on the MK lattice. Similarly, for the WB lattice we 
find 1/v = 0.69 ± 0.35 for the bimodal distribution and 
l/v = 0.68 ± 0.31 for the gaussian distribution. 

In Tab. lIIIBl we summarize the estimates of the critical 
exponents. Once the exponents z, yh and v are known, 
the exponents a and j3 are obtained from the scaling 
relations^ 

a = 2 - (d-z)v, (14) 
P = {d-y h )v. (15) 

The critical exponents obtained on MK lattice for the 
bimodal case are compatible with those obtained by Cao 
and MachtaP^ The exponents yh and z depend strongly 
on the fractal dimension of the lattice, and their value 
for the d = 3 MK lattice are in good agreement with the 
results for the three dimensional Bravais lattice, whilst 
for the d ~ 3.585 WB lattice they are closer to that found 
for the four dimensional Bravais lattice, cf. Tab. |IIIB| 

The value of the exponent v is larger for the WB hier- 
archical lattice than for the MK lattice, as in agreement 
with the behavior in hyper-cubic lattices of increasing 
dimensions. Although the errors bars for this exponent 
are large, we stress that in this case only the WB lattice 
gives a numerical estimate compatible with the result for 
the cubic lattice. 

We, eventually, notice that the Gaussian and bimodal 
exponents are always compatible with each other and 
belong to the same universality class. 

In the RFIM the bond are at first all equal, and the 
disorder is on site. The major reason of moving from 
MK lattices to more complex HL is the need of a better 
treatment of the bond structure. So in the next sections 
we move to models in which the disorder is on the bond 
from the outset. 
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Lattice 


z 


Vh 


1/u 


bimod. MK 


1.491(3) 


2.991(1) 


0.45(23) 


Gauss MK 


1.486(3) 


2.990(1) 


0.44(20) 


bimod. WB 


1.788(2) 


3.575(1) 


0.69(35) 


Gauss WB 


1.787(2) 


3.576(1) 


0.68(31) 


3d MCP 


1.49(3) 


2.988(4) 


0.73(5) 


4d MCP 


1.779(4) 


3.827(1) 


1.280(2) 



TABLE I: Critical exponents at critical fixed distribution for 
RFIM on necklace MK lattice in Fig. [2]and WB lattice in Fig. 
ID The last two rows relate to simulations on Bravais 3D and 
4D hypercubic lattices with bimodal distributed interactions. 



C. Random Bond Ising model in d < 2.5 

In this Section we consider the Ising model with bi- 
modal ± J bond distribution on hierachical lattices mim- 
icking the topology of the square lattice. 

The initial probability distribution for interactions is 

P(Jij,hij, h \i) = \pS(Jii - J) + (1 - PWa + J )\ 



Hhij)s(hl 



(16) 

where J > and p £ [0, 1] gives the probability of a 
ferromagnetic bond. 

For low enough p, this model on regular lattice has 
an antiferromagnetic phase. On hierachical lattices the 
antiferromagnetic order is preserved under PSRG only 
when the rescaling factor b is odd, so that a symmetric 
phase diagram in p O (1 — p) is obtained, where the 
ferromagnetic phase is replaced by the antiferromagnetic 
phase and viceversa. 

To capture some features of Bravais lattices, we shall 
focus on hierarchical lattices with elementary cell more 
complex than the ID-like MK cells. In particular we 
shall consider the cell proposed by NobreJ^S Fig. [5] with 
rescaling factor b — 3 and its extension to b = 5, Fig. [6] 

Nishimori conjecture. Even though the Random 
Bond Ising model in 2D does not display any spin-glass 
phase, this model on hierarchical lattices is an excellent 
play ground to test Nishimori's conjecture! 54 * 55 ^ 

The idea behind the conjecture stems from noting that 
the partition function Z of non-random Ising models on 
self-dual lattices is itself self-dual.^ That is, if expressed 
in terms of the bond Boltzmann factors u±i(J) = e ±J 
(note the use of the reduced parameters), Z is invariant 
under the exchange u±i(J) f-> u± 1 (J), where the m5 =1 (J), 
Fourier transforms of u±i(J), are the dual Boltzmann 
factors As a consequence the critical point of a self- 
dual model is obtained, if unique, by the fix point condi- 
tion u±i(J) = u* ±1 (J), which yields J c = |log(-\/2 + 1). 

In the case of bimodal random bonds the relevant vari- 
ables, within the replica method approach, are the aver- 
aged bond Boltzmann factors Xk(p, J), which correspond 



to the configuration with the spin connected by the bond 
equal to +1 in n — k replicas and —1 in remaining k 
replicas. 

Self-duality is expressed by the invariance of Z n = Z n 
under the simultaneous exchange Xk (j>, J) <-> x* k (p, J) for 
all k. Unlike the non-random case, it is not possible to 
identify the critical point from the fix point condition 
of the duality relations, because the relations with k — 
0, . . . , n a re not satisfied simultaneously. The Nishimori's 
conjectur e! 58 ! 59 ! identifies the multicritical point (pn, Jn) 
with fix point condition for the leading k = Boltzmann 
factor 



Xo(pn, Jn) = Xq(p n , J n ) 



(17) 



on the Nishimori line e~ 2J = (1 — p) /jjjSilfiJD where en- 
hanced symmetry simplifies the system properties signif- 
icantly. This conjecture is proven exact for n = 1, 2 and 
00. In the limit n — > the condition xq = Xq on the 
Nishimori line becomes 



H(pn) = -pN^og 2 {p N ) - (1 - p N )log 2 (l - p N ) 



1 

2 ' 

where the function H(p^) is the binary entropy. 

This c onjecture has proved wrong f or som e systems on 
HL }24i6i| Qhzeki, Nishimori and Berkeil 2 ^! noted that for 
HL a systematic approximation for the multicritical point 
can be obtained by imposing Eq. (171 at each PSRG 
transformations. 

Note that in the two-dimensional case, even though the 
SG phase is absent, the Nishimori point is expected to 
coincide with a critical point, unstable along the phase 
boundary; so, as usual, in the following we identify the 
Nishimori point as the intersection between the Nishimori 
and the critical line. 

Here we test of the original Nishimori conjecture on 
a Folded Square (FS) hierarchical lattice constructed 
through a bond-moving procedure that, at difference 
with the MK bond-moving prescription, retains the lo- 
cal correlation of the bonds, see Figs. [5]and[6j 



1. Numerical results 

In Fig. [8] we show the p,T(— I /J) phase diagram 
obtained from the PSRG analysis in two dimensions on 
the lattice shown in Fig. [6] with rescaling factor b = 5. 

For p — 1 we find the critical temperature is T c = 
2.269185(1), in agreement with the Onsager solution 
T c = 2/log(l + V2) ~ 2.269185. This result, also found 
for the folded square lattice of Fig. [| with b = 3pH an d 
WB with d = 2.32pSl follows from the duality properties 
of the unit cells P 

Nishimori point. The position of the multi-critical 
point expected from Nishimori's conjecture, Eq. ( [l8| ), is 
p N ~ 0.889972. For the cell of Fig. [5] the estimated 
value is p N = 0.8903(2) (2H(p N ) = 0.998(1)^3 w j t h 
T N = 0.9557(18)P 



For the cell of Fig. [6j performing N$ = 
PSRG calculations with a poll of M = 
configurations, we estimate (cf. Tab. [TTJ 



: 20 independent 
10 5 initial bond 



2.5 



PN 

2H( PN ) 



0.8902(1), 
0.9985(1); 



T N = 0.9571(1) 



thus we conclude that the conjecture fails also on this 
more complex hierarchical lattice. 

The folded square cell estimates with b = 3 or b = 5 
turn out to be in good agreement each other and with 
the estimate pn = 0.8905(5) given by the transfer ma- 
trix approach,^ but are slightly larger than those from 
high temperature series expansion, p^ = 0.886(3)p3 and 
Monte Carlo simulation, p N = 0.8872(8)P 

In Tab. [Il] we report the results on critical points for 
an easier comparison. 

Critical slope. An interesting quantity to study 

is the slope of the critical line close to p = lpS 



1 dT c (p) 



T c (l) dp 



(19) 



The Domany's perturbative approach^! yields s — 
2v / 2/[ln(-\/2 + l)] ~ 3.209 for the Ising with ±J bond dis- 
tribution on square lattice. The approach assumes weak 
disorder, i.e., a qualitatively irrelevant disorder that does 
not undermine the existence of the ferromagnetic phase 
at low T and does not change the universality class of 
the PM/FM transition. By computing s it is, thus, ar- 
gued that one can discriminate whether quenched disor- 
der is a relevant perturbation, causing a change in the 
universality classP^ Ohzeki and collaborators^ suggest 
that Domany's method can be applied to any self-dual 
lattice, and then one can probe the relevance of disorder 
from the slope s also for the HL of Figs. [5]and|6] 

The duality approach gives s = 3.27866..=^Tfor the 
lattice in Fig. [5] with 6 = 3. In the b = 5 case of Fig. 
|6j from a best fit of the points close to p = 1 with a 
pool of size M = 5 • 10 6 with N$ — 20 samples we find 
s = 3.316(11). Both values are incompatible with the 
Domany's prediction s ~ 3.209P2 

According to the argument of Ref. [ST] this would im- 
ply that the disorder should be relevant in these cases. 
However, if it is true that a slope equal to Domany's slope 
would be consistent with the hypothesis of irrelevant dis- 
order, the fact that the slope is different from Domany's 
value does not, actually, imply the opposite (i.e., the ex- 
istence of a new "strong disorder" fixed point). Indeed, 
from a direct inspection, no other fix points are found 
because of the introduction of disorder and consequently 
the universality class does not change. Any coupling dis- 
tribution on the transition line tends to become more 
and more peaked under PSRG transformations, tending 
to the pure FM-PM critical fixed point, as shown in Fig. 
|9j Disorder in the two dimensional Ising model does not 
appear to play any role on any lattice cell analyzed here 
(see also the Harris criterion below). 
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FIG. 8: Phase diagram of Ising 2D with "folded square" HL 
with 6 = 5. The diagram is symmetric in p — > 1 — p and 
the anti-ferromagnetic part is not shown. The dashed line 
rappresents the Nishimori line. 




FIG. 9: PSRG evolution of the coupling distribution on the 
critical line near the pure fixed point (steps from 4 to 11 are 
shown for p = 0.01): the width quickly decreases, so that the 
distribution tends to a Dirac delta and the disorder disap- 
pears. No other fixed ditribution. besides the pure criticality, 
are obtained on the PM-FM transition line. 



HL / Method 


Fig./Ref. 


Pn 


FS 6 = 3 


||/2g 


0.8903(2) 


FS 6 = 5 




0.8902(1) 




m 


0.8905(5) 


Square 




0.886(3) 




M 


0.8872(8) 


Nishimori conjecture 


m 


0.889972... 



TABLE II: Multicritical point as computed on different HLs 
and different estimates on the square lattice. 
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Harris criterion. The widely accepted form of the 
Harris criterion^ is that in ferromagnetic systems with 
random interactions the randomness is irrelevant if a, the 
specific heat exponent of the corresponding pure system, 
is negative, while for systems with positive a the random 
system exhibits different critical behavior in presence of 
disorder. For the folded square lattices for both 6 = 3 
and 6 = 5 we find a < 0. 

Though it is known that this criterion can fail on 
HL if the bon ds in the rescaling volume are not all 
equivalent l 44 l 68 tf ZH i n the present case the Harris crite- 
rion confirms that the disorder is, actually, irrelevant, as 
discussed above with respect to the critical slope value. 

FM line reentrance. An important feature of the 
p, T diagram, according to the duality requirements, is 
the reentrance of the transition line below the multicrit- 
ical point: p^ > Pt=o- The zero temperature tran- 
sition point pt=o can be estimated by finite size scal- 
ing analyses of the ground state. Calling E p and E a 
the ground-state energies with, respectively, periodic and 
anti-periodic boundary conditions in one direction, and 
A = E p — E a the domain wall energy, we can determine 
two estimates of critical concentrations of antiferromag- 
netic bonds by looking at the point where the asymptotic 
dependences [A] and [A 2 ] 1 / 2 change from increasing to 
decreasing, where the average [. . .] is taken on different 
bond samples. More explicitly, by defining 



[A] - L p and [A 



2i 1/2 



L 



via the determination of exact ground states for large 
system sizes and huge sample numbers, Kawashima and 



Rieger give the estimates p T 



(i) 



0.896(1) 



and Pyl, 



0.894(2) looking at the point where p and 6 change signf^l 
The PSRG approach with the folded square of Fig. [5] 
(6 = 3) leads to p T =o = 0.8951(3),^ while with the cell 
of Fig. |6] (b = 5) we find p T=0 = 0.8966(2). 

The critical indices for the pure ferromagnetic point 
can be evaluated as discussed in Sec. MI At we obtain 
y T = 0.7303(1) and y h = 0.8518(1) for the folded square 
lattice with 6 = 3, while for the folded square lattice with 
6 = 5 we find y T = 0.7589(1) and y h = 1.059(1). In both 
cases they are not consistent with the values of Onsager 
solution, which read yx = 1 and yh = 1.875, respectively. 
From the knowledge of yx and yh the critical indexes are 
obtained from the usual scaling relations: 



V = 


1 

Vt' 


(20) 


V = 


d + 2-2y h 


(21) 


a = 


2--, 
Vt 


(22) 


P = 


d - Vh 
Vt 


(23) 


7 = 


tyh - d 
Vt 


(24) 


S = 


Vh 
d-yh 


(25) 





FS b=3 


FS b=5 


Exact 


a 


-0.7385(1) 


-0.6353(1) 





P 


1.572(1) 


1.240(1) 


0.125 


7 


-0.4057(1) 


0.1558(1) 


1.75 


6 


0.7419(1) 


1.126(1) 


15 


V 


1.369(1) 


1.318(1) 


1 


n 


2.296(1) 


1.882(1) 


0.25 



TABLE III: Critical indices of pure ferromagnetic critical 
point for Ising model on the folded square lattice for 6 = 3, 
Fig. [5] and for 6 = 5, Fig. [6] 



Lattice type 


F 


k 


Ref. 


e 


MK b — 2 




M 


-0.270(2) 


MK b = 3 




[76] 


-0.278(2) 


WB b = 2 




[34] 


-0.290(3) 


WB b = 3 






[76] 


-0.298(2) 


FS 6 = 3 




5 


133 


-0.275(1) 


FS 6 = 5 




6 




-0.2714(2) 






m 


-0.291(2) 


Square 




m 


-0.287(4) 








-0.284(4) 



TABLE IV: Stiffness exponent 6 on different HLs and different 
estimates on the square lattice. 



and their numerical values are reported and compared 
on Tab. 



Ill 



(141, 



Notice that, at difference with Eq. 
here the fixed point is at finite T (as 1/J) and the index 

z = oP 

Zero temperature stiffness. We conclude this sec- 
tion by discussing the exponent v of the zero-temperature 
spin-glass transition for the case of a Gaussian distribu- 
tion of bonds with zero mean and initial width oj. It can 
be obtained directly from scaling of ctj under PSRG : 



,(b) ~ ajb 6 



(26) 



The sign of the stiffness exponent is directly related 
to the low temperature phase: for positive (negative) 
9 the system scales under PSRG flow towards strong 
(weak) couplings, distinctive of a low temperature spin- 
glass (paramagnetic) phase. For continuous and sym- 
metric probability distributions P(J), the temperature T 
appears in the PSRG equations as a dimensionless ratio 
between couplings, so that the scaling (261 is equivalent 
to T - L e , or L ~ T 1 / 6 . In a phase transition at T ->• 
the latter scaling can be identified with the scaling of the 
correlation length £ ~ T~ v implying!^ 



v=-\i 



(27) 



In Fig. [TO] it is shown the behavior of aj as function 
of PSRG step for the case of zero-average Gaussian ini- 
tial bond distribution on the folded square cell of Fig. 
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5 10 15 20 25 30 

RG step 

FIG. 10: Standard deviation of the Gaussian probability dis- 
tribution of quenched disorder at T = vs. PSRG iteration 
steps. 




FIG. 11: Cubic cells of length 6 = 3 on the left hand side, 
and relative folded cube hierarchical lattice on the right hand 
side. The two roots are open circles. All incoming sites and 
all outcoming ones, pointed by the arrows in the left figure, 
are put together and generate root sites 1 and 2 on the right 
hand side figure. The inner sites of the cubic cell become the 
inner sites of the hierarchical cell. 



[6j From this we get the 6 = -0.2714(2), leading to 
v = 3.685(3). For the cell of Fig. [5l with b = 3 it was 



found 9 = -0.275(1) P We report in Tab. [fv]the val- 
ues obtained on different HL d 25 * 34 * 46 * 76 ! and on the square 
Bravais lattice.^^U We observe that the value for the 
folded square with both b = 3 or b = 5 is similar to that 
found for the MK lattice, whilst for the WB lattices the 
values are closer to that of the regular lattice, especially 
the case with b = 2. 

Passing from b = 3 to b = 5, and thus increasing the 
connectivity of the lattice, we obtain an improvement, as 
compared to square lattice, only for the pure model (al- 
though the estimates for the critical exponents are still 
not compatible with the exact values of the square lat- 
tice), whilst the values for the disordered quantities are 
generally compatible with each other. 

To obtain tests for quantities with relevant disorder 
we have to move to models in which there is a spin-glass 
phase, which we will analyze in the following sections. 



D. Random bond Ising model in d > 3 

In this section we consider the Ising model with the ± J 
coupling distribution, Eq. ( 16 ) on hierarchical lattices 
mimicking the topology of the cubic lattice. Recently, 



Salmon, Agostini and Nobre have studied the Ising spin 
glass on the hierarchical Wheatstone bridge (WB) lat- 
tices drawn in Figs, [l] [4^51 obtaining accurate phase 
diagrams and showing that on these lattices the lower 
critical dimension for the spin glass phase is greater than 
d = In 5/ In 2 ~ 2.32. The next pattern of the WB family, 
cf. Fig. [4j corresponding to a lattice in dimension higher 
than 2.5 has a fractal dimension d ~ 3.58. Large devia- 
tions from the cubic lattice behavior might, then, occur 
and results somewhat in the middle between 3D and 4D 
Bravais lattice models. 




61 81 



FIG. 12: MK lattice with b = d = 3. 



We, here, investigate the model on the Folded Cube 
(FC) hierarchical lattice shown in Fig. 



11 This is the 



"three-dimensional" extension of the folded square lattice 
of Fig. [5j and it was introduced to study the anisotropic 
ferromagnetic Potts model in three dimensions.^! This 
lattice has a fractal dimension d = In 35/ In 3 ~ 3.2362, 
always larger than 3, though nearer to it as compared 
to the WB. It has 6 = 3 and, unlike the latter, it is 
able to retain a possible antiferromagnetic order, as well. 
As a further comparison we will report on the critical 
properties of the Ising spin-glass model on a MK lattice 
with b = 3 and fractal dimension d — 3, cf. Fig. [L2| 
introduced in Refs. |35l81j . 

The resulting phase diagram is shown in Fig. |13| An 
important feature of the phase diagram is the small rcen- 
trance in the region below the multicritical point. As a 
consequence, by lowering the temperature one goes from 
the high temperature (disordered) paramagnetic (PM) 
phase to an ordered ferromagnetic (FM) phase and, even- 
tually to a low temperature disordered spin-glass (SG) 
phase. The values of the critical points are reported in 
Tab. El 
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FIG. 13: Phase dia gram of Ising spin glass model on the folded 
cube HL, cf. Fig. [ll] obtained for M = 5 ■ 10 4 e N s = 10. 
The dashed line represents the Nishimori line. The plot is 
symmetric in p — > 1 — p and an antiferromagnetic phase is 
present a small p, in place of the FM one. 



Fixed points (p, T) coordinates 



HL 


FM 


SG 


T = 


MC 


Fold cube 


1 


0.5 


0.764(2) 


0.7547(3) 


Fig. [I] 


5.066(1) 


1.072(1) 





1.779(1) 


WB "3D" 


1 


0.5 


0.760(1) 


0.745(2) 


Fig. |4l 


5.457(1) 


1.112(2) 





1.620(2) 


WB "4D" 




0.5 


0.667(2) 


0.664(2) 


Ref. [25] 




2.515(2) 





2.836(2) 


MK 


1 


0.5 


0.761(1) 


0.752(7) 


Fig. m 


5.383(1) 


1.136(4) 





1.797(3) 



TABL E V : Position of typical critical points of phase diagram 
in Fig. 13 for the hierarchical lattice 11 ferromagnetic (FM), 
spin-glass (SG), zero temperature (T = 0) and multi critical 
point (MC). For comparison also values of the same fixed 
points are reported for other hierarchical lattices discussed in 
the text. 



1 . FM fixed point 

The transition for p = 1 on the folded cube cell is ob- 
tained at T c = 5.066(1), 12% larger as compared to the 
same model on cubic lattice (where T c = 4. 5115. .P). On 
the d ~ 3.58 WB lattice the difference was about 21%P 
On the MK lattices the best known result is obtained for 
the d — 3 lattice in Fig. 12 where the critical temper- 
ature of the pure transition is T c = 5.38(3), 19% larger 
than the one in the cubic latticeP 

In the folded cube the points on the critical line be- 
tween the FM and the PM phases are attracted to a sin- 
gle fix point located at p — 1. The critical exponents at 



critical 
index 


MK 
Fig. [li] 


WB 
Fig. [4] 


Folded cube 
Fig. [ll] 


Cubic 
Ref. El 


a 


-0.2169(1) 


-1.121(1) 


-0.1253(1) 


0.110(1) 


P 


1.112(1) 


2.282(1) 


0.9012(1) 


0.3265(3) 


7 


-0.006490(1) 


-1.443(1) 


0.3229(1) 


1.2372(5) 


5 


0.9942(1) 


0.3676(1) 


1.358(1) 


4.789(2) 


V 


0.6850(1) 


0.8705(1) 


0.6567(1) 


0.6301(4) 


V 


2.009(1) 


3.658(1) 


1.508(1) 


0.0364(5) 



TABLE VI: Critical exponents in pure ferromagnetic critical 
point for Ising model on the folded cubic cell and comparison. 



this point are yr — 1.523(1) and yh = 1.864(1), identify- 
ing a second order transitionPS On other HLs we obtain: 
y T = 1.149(1) and y h = 0.9636(1) for the WB lattice 
in Fig. [4] y T = 1.460(1) and y h = 1.613(1) for the MK 
lattice in Fig. [12J These are to be compared with nu- 
merical estimates by means of simulations on the cubic 
lattice: y T = 1.587(1) and y h = 2.482(1) P Though none 
of them agrees with the cubic lattice results the Folded 
cube lattice in Fig.[Tl]yields the nearest estimate for both 
exponents. The corresponding physical exponents are re- 
ported in Table |VI| Note, in particular, that for all the 
HLs considered here we have a < and the disorder is 
irrelevant, thus satisfying the Harris criteriorP: the dis- 
tributions on the critical PM-FM line flow to the pure 
model critical point. 



2. Spin-glass fixed point 

For the DC HL in the totally disordered case p = 0.5 
the PM-SG transition is located at T c = 1.072(1), with 
a difference of 4.5% relatively to the Bravais cubic lat- 
tice critical point, for which the Monte Carlo simulations 
provide T c = 1.120(4)P As reported in Tab. [vjthe WB 
lattice in Fig. [I]has T c = 1.112(2) and the MK lattice in 
"has T c = 1.136(4). 
oehavior of ln(crj) under successive PSRG itera- 
0.5 is shown in Fig. [Ml From its linear 



Fig. 12 



The 
tions for p 

behavior we estimate for the stiffness exponent of the 
spin-glass phase = 0.2052(1) (cf. Sec. |III C 1[ ). For the 
cubic lattice one finds in the literature = 0.19(1^^ or 
9 = 0.20(5)P We, thus, obtain a result very close to 
that expected for the cubic lattice. On the WB lattice 
one has 9 = 0.297(3)P while for the MK lattice it is 
found 9 ~ 0.27P The result are summarized in Table 

rvm 

All critical points along the PM-SG transition are at- 
tracted by an unstable fixed point at p = 0.5. In Fig. [15] 
we show the unstable fix point coupling distributions for 
the three lattices considered so far. Critical exponents for 
this (unstable) fix distribution are computed by adapt- 
ing the methods used for the RFIM in Sec. |IIIB| A first 
exponent y^ can be obtained generalizing the pure case 
definition: A/j = {dhhn), where the average is carried 
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Fig 


. Ref. 
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MK d= 3 


12 


m 


0.27(1) 


WB d~ 3.58 


4 


m 


0.297(3) 


FS d~ 3.24 


11 




0.2052(1) 


Cubic 


m\ 


0.19(1) 
0.20(5) 



TABLE VII: Stiffness exponent as computed on different HL 
and different estimates on the cubic lattice. 



out by extracting the interactions from the fix point dis- 
tribution. Indeed, in our case this is easy and can be 
exactly calculated. In fact for the fix point distribution 
we have h = h) = 0, and so H(s) = H(—s). It is, then, 
straightforward to obtain, using Eqs. Q and Q that 
dhfiR = c for every choice of the J interactions, where c 
is the number of internal sites connected to each external 
site: c = 9 for the folded cube cell in Fig. [II] and the 
MK in Fig. [l2j c = 4 for the WB in Fig. g] For all these 
cells the exponent y^ is then — log b (dhhji) — 2 (see 



Sec. Ill A). We note that yh < d, ensuring that the mag- 
netization is continuous at the transition (no first order 
transition) . 

To get the correlation length exponent v we note that 
defining t = aj — <jj, where cr} is the value at the critical 
point, we obtain the scaling law (note that, as always in 
this paper, we are using the reduced variables /3 J — > J) 



(28) 



By taking n such that b nx t — to, where to is arbitrary 
but fixed, we end up with 




FIG. 15: Fixed probability distribution for the spin glass - 
paramagnet transition for Ising spin glass on lattice in Fig. 
T2](MK), Fig. [|(WB) and Fig. [ll](FC). These distributions 
are reached after about 2 — 3 steps of renormalization, starting 
from the binomial distribution, and maintain this shape for 
about 10 — 15 steps steps, after which they starting moving 
away because of statistical fluctuations. 



Lattice 


Fig. 


Resc. 


Dim. 


1/1/ 


MK 


12 




b = 3 


d = 3 


0.297 ±0.026 


WB 


4 




b = 2 


d ~ 3.58 


0.308 ±0.063 


Folded cube 


11 




6 = 3 


d ~ 3.24 


0.262 ±0.037 


Bravais cubicp^ 








d = 3 


0.408 ± 0.025 



TABLE VIII: Estimates for the exponent v at transition SG- 
para for the Ising spin glass model. 



£(*) = 



i(to) ~ t- 



(29) 




FIG. 14: Graph of ln(crj) under repeated applications of the 
PSRG transformation at zero temperature and p = 0.5 on 
lattice in Fig. [TT] 



from which we get v = 1/x: the value of v can be esti- 
mated by the trend of aj near the critical point distri- 
bution. 

We use the method already exposed in Sec. |III B| after 
obtaining an instance of the fixed distribution, we make 
a copy of it and multiply each coupling of the copy by 
a quantity (1 + 5) sa 1.05. The two copies of the distri- 
bution are, then, simultaneously iterated in the PSRG 
transformation, with the first copy forced near the fixed 
point and the second one free to flow. We eventually es- 
timate v by means of Eq. (13), where the parameter t n 



is now the difference between the values of <jj in the two 
copies at RG step n. 

Continuing this way we find the values shown in Ta- 



ble VIII for different HLs and we notice that they are 
compatible with each other within the statistical error. 
When compared to the estimate for the cubic Bravais 
lattice, though, only the WB one turns out to be nearly 
consistent with v = 2.45(15) ->• Xjv = 0.408(25) P 

Proceeding as for the correlation length, and using the 
free energy / scaling law 



/(*) = b~ d f(t') = b- d f(b x t) 



-nd 



f{b nx t) 
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critical 
index 


MK 
Fig. [li] 


WB 
Fig.0 


Folded cube 
Fig. [n] 


Cubic 
Ref. E] 


a 


-8.10(88) 


-9.6(24) 


-10.4(17) 


-5.4(5) 




3.37(30) 


4.20(86) 


4.27(60) 


0.77(5) 


7 


3.37(30) 


1.35(28) 


2.92(41) 


5.8(4) 


V 


3.37(30) 


3.25(66) 


3.82(54) 


2.45(15) 


V 


1 


1.585... 


1.236... 


-0.375(10) 



TABLE IX: Estimates for the physical exponents at transition 
SG-para for the Ising spin glass model. 




with n such that b nx t — to and to arbitrary but fixed, we 
have 



FIG. 16: MK lattice with fractal dimension d = 3.58, the 
same as WB lattice in Fig. [3] 



(30) 



Then from d 2 f/dt 2 ~ \t\ Q , we obtain the scaling rela- 
tion 



as well as 



where 



a = 2 = 2-dv, 

x 



= v{\ + ij) 

2 

7 = (2-f?K 



rj = d+2-2y h = d-2 



(31) 



since = 2 for all our HLs. The estimates of the physi- 
cal exponents are reported in Tab. |IX| 

Concluding this section, on the FC lattice in Fig. [TT] 
the estimates for the pure criticality are much closer to 
those on cubic lattice (although still not compatible), 
compared to MK in Fig. [12] and WB in Fig. |4j Also 
for the stiffness exponent of SG phase a remarkable im- 
provement is observed: its estimate is compatible with its 
cubic value on FC lattice but on on other HLs. For others 
quantities, though, at the disordered fixed point such im- 
provement is unseen. In particular in the estimate for the 
v exponent at the SG-PM transition. We can argue that 
this is so because the stiffness exponent depends more 
on the local geometrical properties of the lattice, that is 
really improved moving to the FC lattice, whilst for the 
critical properties depends on longer distances, that are 
still dominated by the hierarchical backbone. 

More stringent tests can be obtained for the Blume- 
Emery- Griffiths model, that we will analyze in the next 
section. For this system the inverse first order transi- 
tion expected from mean-field theory^ and simulation 
in finite dimension^ is absent on the MK lattice in Fig. 

0]p] 



IV. BLUME-EMERY-GRIFFITHS MODEL 

We now move on to a different system, the Blume- 
Emery-Griffiths (BEG) model, originally devised to 
study the superfluidity transition and phase separation 
in He 3 -He 4 mixtures.^ The model is known to display, 
besides a second order phase transition, also a first order 
transition, both in the ordered case (between the PM and 
FM phases) and in the case with quenched disordered in- 
teractions (between the PM and SG phases). 

The ordered model cases have been introduced and 
solved in the mean-field approximation in Refs. |88l - 
190] - Finite dimensional analysis has been carried out by 
different means, e.g., series extrapolation techniques,^ 
PSRG analysis,^ Monte Carlo simulations,^ effective- 
field theory^ or two-particle cluster approximation.^ 

Extensions to quenched disordered models, both per- 
turbing the ordered fixed point and in the regime of 
strong disorder, have been studi ed throu ghout the years 
by means of mean-field theory ffl 95 * 96 ! PSR G analysis 
on Migdal-Kadanoff hierarc hical lattic es 1251451 and Monte 
Carlo numerical simulations! 31 * 37 * 38 * 97 ! The latter studies 
show that a critical transition line separates the SG and 
PM phases. Like in the mean-field cases, it consists of a 
second order transition terminating in a tricritical point 
from which a first order inverse transition starts! 3 ^ Fur- 
thermore, a reentrance of the first order transition line is 
present for positive, finite values of the chemical poten- 
tial of the holes' yielding the so called inverse freezing 
phenomenon of an amorphous phase arresting itself in a 
blocked, solid-like state upon heating. 

In absence of any purely ferromagnetic contribution, 
the PSRG approach on MK cells apparently does not 
show any first order phase transitions, nor any reen- 
trance, as shown by Ozcelik and BerkerP^In this section 
we will deepen their analysis at p = 1/2, investigating 
the critical behavior on the MK lattices of Fig. [12] and 
Fig. 16 on the Wheatstone Bridge (WB) lattice of Fig. 
[4] and on the Folded Cube (FC) lattice of Fig. 



11 



The BEG model with generic magnetic exchange in- 
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teraction is defined by the Hamiltonian (we use reduced 
variables) 



mis) = J2 + K J2 ^ - A E 



(32) 



<tj> 



<tj> 



where Sj = ±1,0 and Jy can be deterministic or 
quenched disordered and distributed according to some 
probability distribution. In the latter case, under 
PSRG transformation all renormalized interactions be- 
come quenched random and it is convenient to start the 
iteration using the more general form 



- 2 s 2 

■ 3 



<*i> («') 



(33) 



The model is, further, defined, by the multivariable initial 
probability distribution of the interactions: 



(34) 



x«5(i^) 5(Atj - A ) <J(Aj,) 



We notice that if Ao <C —1 the values s, = are sup- 
pressed and the model tends to the Ising spin glass model 
analyzed in previous sections: 

-m{{Jij}; {>*}) = 53 Jij s i s 3 

with Si = ±1 and 

P( J -) = J dK tJ d\, dA\ 3 VU,, , K tj , Ay , A^.) . 

Decimating the inner sites at a given hierarchical cell 
with fixed outer sites s a and Sj,, and using the up-down 
symmetry of the Hamiltonian, the relations for the renor- 
malized interactions imposed by the conservation of the 
partition function, cf. Eq. Q, can be written similarly 
to Eqs. (|5j>-(j4j> as 



Jr 
Kb. 
Ah 



1 , ( x++ 

2 log U+r 



hog( x+ \ x+ - x ° o 

2 V x o+ x +o 



1 00 



2 \x +0 x 0+ 



1 , / x +o 

-log 

2 V x + 



(35) 



where x SaSb are the Boltzmann factors, cf. Eq. Q. In 




FIG. 17: Phase diagram for the BEG model on the MK(d = 3) 
lattice in Fig. [TJl WB lattice in Fig. [4] FC lattice in Fig. [TT] 
and simulations . The line for the MKs and WB are obtained 
with M = 10 6 , whilst for the FC we use M — 5 ■ 10 4 . Inset: 
detail of the reentrance region for WB and 3D cubic lattices, 
compared to the MK HL line (no reentrance). 



the T — limit these relations become 



2J R 
2K R 



max[-/3-H(l, 1, s)] - max[-^(l, -1, s)] 
max[-/3ft(l, 1, s))] + max[-m(h -1, »)] 
- 2{ max[-/9H(l, 0, s)] - max[-/3ft(0, 1, s)] } 
+ 2max[-/3"H(0,0,s)] (36) 
2A fl = 2max[-/3ft(0,0,s)] 

- max[-/3H(l,0,s)] - max[-/3H(0, 1, s)] 

max[-mO-, 0, s)] - max[-0H(0, 1, «)] 



2A 



The phase diagrams relative to the different hierarchi- 
cal lattices are shown in Fig. [17] They are obtained 
representing the probability distributions by a pool of 
M = 10 5 interaction quadruples ( J, K, A, At) and every 
evolution is repeated N$ — 10 times. As worked out 
in Sec. Ill paramagnetic, ferromagnetic and spin glass 



phases are determined from the analysis of the PSRG 
flux of J = (Jn) and 



„ = <4>-<^> 2 - 

In all cases we obtain a second order transition between 
paramagnetic and spin glass phases, with all points on 
the transition line attracted by a unique fixed distribu- 
tion at A -> — oo (see Fig. 18), thus belonging to the 
same universality class of the SG-PM transition in the 
Ising spin glass investigated in the previous section. 

Also studying hierarchical lattices much more complex 
than MK, the first order transition typical of the BEG 
model on Bravais lattice is not found, so we have a strong 
indication that this is an intrinsic limit of hierarchical 
lattices, and not only of the MK kind of cells. 

It is important to emphasize, however, that on the WB 
d w 3.58 lattice we obtain a feature not present for FC 
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FIG. 18: (Projection of unstable fixed distribution 
V{ J, K, A, A*) for BEG model at the transition SG-para on 
the WB lattice in Fig. [1] In the distribution A is runaway 
to minus infinity exponentially, while the other interactions 
remain finite and the distribution retains this shape. The 
shape of relative distribution on the MK and FC lattices is 
very similar. 



lattice and for MK lattices, neither in fractal dimension 

3.58 (cf. Fig. 



3 (cf. Fig. 12 1 nor in d ~ 3.58 (cf. Fig. |16|: an inverse 
transition between spin glass and paramagnet. Evidence 
for the reentrance is clearly obtained using the T = 
Eqs. (|36|). Inverse freezing is predicted in mean-field 



theory 3 ^ and found in 3D numerical simulations.^ 



V. CONCLUSIONS 

We have provided a review of the standard methods 
to develop PSRG on hierachical lattices. We stress how 
on one side this analysis can be useful to test general re- 
sults (e.g. Nisfiimori conjecture, Harris criterion), and on 
the other side it is far easier and faster than Montecarlo 
simulations on Bravais lattice. 

We have used these methods to investigate the Ran- 
dom Field Ising Model (RFIM) , the Random Bond Ising 
Model (RBIM) and the Blume-Emery-Griffiths (BEG) 
model defined on several MK and non-MK hierachical 
lattices, obtaining the whole phase diagrams and critical 
exponents. 

The RFIM has been analyzed on a non-MK (WB in 
Fig. [4j for the first time, and we have shown that the bi- 
modal and Gaussian disordered cases belong to the same 
universality class. 

The RBIM in d < 2.5 has been analyzed on the folded 
square HLs family for 6 = 3 (Fig. [5]) and 6 = 5 (Fig. [6]), 
where we find that the Nishimori conjecture fails and the 
disorder is irrelevant. 

The RBIM in d > 2.5 has been analyzed on the folded 
cube (Fig. Ill and compared to MK (Fig. 



SG-PM transition and also show that the Harris criterion 
holds. 

The BEG model has been analyzed on non-MK HLs 
for the first time, WB (Fig. [§ and folded cube (Fig. [ll]), 
and compared to MK with d — 3 (Fig. 12) and d ~ 3.58 
(Fig. 16), and in al l the cases the first order transition 
of the Bravais lattice^^H is absent. 

Our results provide a clue to the possibility of obtain- 
ing approximations of models on regular lattice by a the 
similar models on hierachical lattice. 

We show that it is possible to obtain a good picture 
of the actual phase diagram, but more difficult a proper 
determination of critical exponents. 

By introducing more complex elementary cells, with 
some non trivial internal structure, one hopes of captur- 
ing the local geometrical properties of the bonds. At least 
part of it. For pure models, although it is not a system- 
atic approximation, a general improvement is obtained 
using more complex unit cells, that locally mimic better 
the connectivity of the Bravais lattice. In the disordered 
case, instead, the situation is less definite, and no net 
improvement is observed: the WB (Fig. |4| proves to be 
the slightly most reliable (generally quantitatively better 
than the more complex folded cube in Fig. 11), and in 



particular shows the expected inverse transition for the 
BEG model, but we can not give a general explanation 
for this. 

In particular, the fractal dimension seems to play a 
minor role, as the three dimensional regular lattice is 
better approximated by the WB with d ~ 3.58, compared 
to d = 3 MK (Fig. [l1| and d ~ 3.24 folded cube (Fig. 



Ill, while the d ~ 3.58 MK in Fig. 16 is the worst 



12j) and WB 

(Fig. [4j, where we obtain the critical exponents for the 



approximation by far. 

On the other hand, the scaling factor 6 has the known 
role to determine if the antiferromagnetic order can be 
preserved, as it is possibile only when 6 is odd so that 
negative interactions in the unit cell involve negative in- 
teractions between the external sites (and this leads to 
a phase diagram symmetric under the inversion of the 
bonds sign). Our results indicates that this feature does 
not play a crucial role in the disordered systems (at least 
until the negative bonds become dominant). The two in- 
vestigated HLs with an even 6, the WB and the d ~ 3.58 
MK, indeed, appear to be, respectively, the best and the 
worst in approximating the phase diagram of the models 
on regular lattice. The only case in which the folded cube 
lattice provides a remarkable improvement with respect 
to the WB lattice is in the estimate of the SG stiffness 
exponent. A more structured inner connectivity could, 
thus, become important at low temperatures. 

In conclusion, our results show that the approximation 
on HL is particularly poor for disordered systems, and 
strongly suggest that its limitations are intrinsic to the 
hierarchical nature. The most striking case is the lack of 
first order transition in BEG. Indeed, using more complex 
HLs we may have a better treatment of short distances, 
i.e., short loops, but longer distances seem to be definitely 
dominated by the hierarchical backbone. 
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